Global longitudinal strain and plasma biomarkers for prognosis in heart failure complicated by diabetes: a prospective observational study

Background Heart failure (HF) and diabetes are associated with increased incidence and worse prognosis of each other. The prognostic value of global longitudinal strain (GLS) measured by cardiovascular magnetic resonance (CMR) has not been established in HF patients with diabetes. Methods In this prospective, observational study, consecutive patients (n = 315) with HF underwent CMR at 3T, including GLS, late gadolinium enhancement (LGE), native T1, and extracellular volume fraction (ECV) mapping. Plasma biomarker concentrations were measured including: N-terminal pro B-type natriuretic peptide(NT-proBNP), high-sensitivity troponin T(hs-TnT), growth differentiation factor 15(GDF-15), soluble ST2(sST2), and galectin 3(Gal-3). The primary outcome was a composite of all-cause mortality or HF hospitalisation. Results Compared to those without diabetes (n = 156), the diabetes group (n = 159) had a higher LGE prevalence (76 vs. 60%, p < 0.05), higher T1 (1285±42 vs. 1269±42ms, p < 0.001), and higher ECV (30.5±3.5 vs. 28.8±4.1%, p < 0.001). The diabetes group had higher NT-pro-BNP, hs-TnT, GDF-15, sST2, and Gal-3. Diabetes conferred worse prognosis (hazard ratio (HR) 2.33 [95% confidence interval (CI) 1.43–3.79], p < 0.001). In multivariable Cox regression analysis including clinical markers and plasma biomarkers, sST2 alone remained independently associated with the primary outcome (HR per 1 ng/mL 1.04 [95% CI 1.02–1.07], p = 0.001). In multivariable Cox regression models in the diabetes group, both GLS and sST2 remained prognostic (GLS: HR 1.12 [95% CI 1.03–1.21], p = 0.01; sST2: HR per 1 ng/mL 1.03 [95% CI 1.00-1.06], p = 0.02). Conclusions Compared to HF patients without diabetes, those with diabetes have worse plasma and CMR markers of fibrosis and a more adverse prognosis. GLS by CMR is a powerful and independent prognostic marker in HF patients with diabetes. Supplementary Information The online version contains supplementary material available at 10.1186/s12872-024-03810-5.


Background
The Asian continent has the highest prevalence of heart failure (HF) cases globally [1].Diabetes is especially common in South East Asian HF populations, where a unique 'lean-diabetic' phenotype with worse outcomes has been identified [2].The Asian Sudden Cardiac Death in Heart Failure (ASIAN-HF) registry (across 11 Asian countries) reported a 42.5% prevalence of diabetes among HF patients, specifically in higher-income countries such as Singapore, Hong Kong, and Thailand [3].Patients with HF and diabetes have consistently worse clinical outcomes, including higher risk of hospitalisation for HF and death, compared to those without diabetes [4][5][6].These findings appear to hold regardless of whether the HF is ischaemic or non-ischaemic in etiology, and regardless of left ventricular (LV) ejection fraction [5,6].
Cardiac magnetic resonance (CMR) imaging has become the non-invasive reference standard for evaluating HF due to its ability to accurately assess cardiac morphology, function, and myocardial tissue characteristics.In particular, late gadolinium enhancement (LGE) permits visualisation of focal replacement myocardial fibrosis, while T1 mapping pre-and post-gadolinium contrast enables non-invasive measurement of the myocardial extracellular volume fraction (ECV), a quantitative marker of myocardial diffuse interstitial fibrosis.Additionally, global longitudinal strain (GLS), defined as the change in the LV myocardial length between diastole and systole divided by the original end-diastolic length, provides a measure of LV systolic function by CMR that is effectively the same as GLS measured by echocardiography [7,8].
Patients with diabetes frequently have impaired GLS and have a higher degree of myocardial fibrosis as assessed by ECV and histology [9,10].GLS by CMR appears to have prognostic value in HF regardless of ejection fraction and whether the cause of HF is ischemic or non-ischemic [11][12][13].There is a growing body of literature demonstrating the prognostic value of GLS in asymptomatic patients with diabetes [14,15].However, the prognostic utility of GLS by CMR in diabetes patients with established HF is unknown.
This study aimed to assess the prognostic significance of GLS by CMR and novel HF plasma biomarkers associated with inflammation and fibrosis in a cohort of patients with HF and diabetes across the spectrum of LV ejection fraction.We hypothesized that GLS would have an incremental prognostic association in this group, beyond plasma HF biomarkers, LV ejection fraction, and CMR markers of myocardial fibrosis.

Study population
Patients with HF were recruited prospectively across six tertiary cardiac centres in Singapore (Asian neTwork for Translational Research and Cardiovascular Trials [ATTRaCT], ClinicalTrials.govNCT02791009).Patients were included if they presented to hospital with a primary diagnosis of HF, or if they attended a hospital clinic within 6 months of an episode of decompensated HF (requiring hospitalization or treatment in an out-patient setting).In all cases, a trained cardiologist adjudicated the clinical diagnosis of HF.The exclusion criteria were: HF primarily due to severe valve disease, HF due to acute coronary syndrome resulting in a transient episode of acute pulmonary oedema, severe renal failure (estimated glomerular filtration rate < 15 mL/min per 1.73m 2 ), specific causes of HF (constrictive pericarditis, complex adult congenital heart disease, hypertrophic cardiomyopathy, eosinophilic myocarditis, cardiac amyloidosis, and acute chemotherapy-induced cardiomyopathy), isolated right HF, and life threatening non-cardiac co-morbidity with life expectancy of < 1 year.All patients underwent clinical assessment at baseline.Diabetes status was identified by baseline questionnaire at recruitment.

Biomarkers
Blood was collected in dipotassium (K 2 )-EDTA vacutainer tubes and transported on ice for processing within one hour.Plasma was separated by centrifugation at 3500 g for 10 min at 4 o C and stored at -80 o C until analysis.

Clinical outcomes
The primary outcome was a composite of time to either first hospitalization for HF (regardless of prior history of hospitalization for HF) or all-cause mortality.Follow-up was conducted through a clinic visit at 6 months after baseline assessment and structured phone interviews with the participants at the 1-and 2-year timepoints.Data in patients who were lost to follow-up were censored at the date when the patient was last known to be alive and had not experienced an event.

Cardiovascular magnetic resonance image acquisition
All patients in the ATTRaCT cohort were invited and assessed for suitability for CMR.Those who agreed and were eligible underwent a standardized CMR protocol with a 3 Tesla MRI scanner (Ingenia, Philips Healthcare, Best, The Netherlands).Balanced steady-state free precession cines were acquired in the standard longaxis views and a short-axis stack from base to apex, as described previously [20].LGE images were acquired at 10 min after 0.1 mmol/kg of gadobutrol (Gadovist®, Bayer Pharma AG, Germany) with a phase sensitive inversion-recovery gradient-echo imaging sequence.Typical parameters were: repetition time (TR) = 6.1 ms; echo time (TE) = 3 ms; time to inversion (TI) = 320-340 ms, flip angle 25°, voxel size = 1.5 × 1.7 × 8 mm, SENSE factor = 2.4, slice thickness 8 mm with 2 mm gap to match shortaxis cine slice positions.The inversion time for optimal myocardial nulling was selected from an inversion time scout sequence.T1 maps were acquired at the basal and mid-ventricular short-axis levels, pre-and 15-min postcontrast with modified Look-Locker Inversion-recovery (MOLLI) 5s(3s)3 and 4s(1s)3s(1s)2s acquisition schemes, respectively [21].

CMR analysis
Image analysis was performed using CVI42 software (Circle Cardiovascular Imaging, Calgary, Canada) by trained imaging fellows at the National Heart Research Institute of Singapore CMR Core Laboratory, who were blinded to the clinical information of the patients.Ventricular volumes, mass, and ejection fraction were measured from the short-axis cine stack, using manual contouring of the left ventricle in end-diastole and endsystole, excluding papillary muscles, as detailed previously [20].LV volumes and mass data were indexed to body surface area.The presence of LGE was assessed qualitatively by two readers according to the recommendations by the Society of CMR [22].Average native and post-contrast myocardial T1 values were measured by placing a region of interest (ROI) within the middle third of the short-axis myocardial wall at the basal-and mid-ventricular levels, while avoiding regions of focal LGE.The myocardium-blood pool interface was carefully avoided in order to prevent partial volume effects.Pre-and post-contrast blood T1 values were measured in a ROI drawn within the blood pool.ECV was calculated from the pre-and post-contrast average blood and myocardial T1 values, as described previously [23,24].Myocardial strain was analysed in the cine images using the Tissue Tracking Plugin [25].

Statistics
Normality was assessed for continuous variables using the Shapiro-Wilk test.Normally distributed data are presented as mean ± standard deviation.Non-normally distributed data are presented as median [interquartile range].Comparisons were performed for continuous variables using the parametric Student t-test or the nonparametric Mann-Whitney U test.Categorical variables are presented as number (percentage) and compared using the χ 2 test.
Univariable Cox regression analysis was performed to identify prognostic variables in the entire cohort and also in the diabetes subgroup.Clinically relevant variables (age, sex, body mass index, diabetes, smoking, coronary artery disease, hypertension, New York Heart Association functional class, and systolic blood pressure), circulating biomarkers (NT-proBNP, hs-TnT, GDF-15, sST2, and Gal-3) and CMR markers of function, remodelling and fibrosis (LV end diastolic volume index, mass index, ejection fraction, GLS, LGE, Native T1, and ECV) were tested in the Cox models.Covariates with a p-value < 0.05 in univariable analyses were entered into the multivariable Cox model to identify independently prognostic variables, using forward stepwise selection (probability to remove threshold ≥0.1).Event-free survival curves were examined using the Kaplan-Meier method and compared with the log-rank test.In the Kaplan-Meier analysis, GLS was treated as a binary variable according to the cohort median, in order to maximise statistical power.Statistical analyses were performed using SPSS Version 28 (Statistical Package for the Social Sciences, International Business Machines, Inc., Armonk, New York, USA) and GraphPad Prism 9.4.1 (GraphPad Software, Inc., San Diego, California, USA).A two-sided p-value < 0.05 was considered as statistically significant.

Results
Figure 1 shows a flow chart describing patient inclusion.Of the 623 patients enrolled in the ATTRaCT study, 523 patients underwent a baseline CMR study.We excluded 168 studies performed without contrast.An additional 40 patients were excluded due to incomplete baseline clinical or CMR data sets, alternative diagnosis on the basis of CMR, or loss to follow-up.In total, 315 subjects (diabetes, n = 159; without diabetes, n = 156) were included in the study cohort.
The baseline clinical characteristics are shown in Table 1.Compared to patients without diabetes, patients with diabetes were older (60±10 vs. 56±12 years, p < 0.001), more likely to have a history of hypertension (76 vs. 52%, p < 0.001), coronary artery disease (72 vs. 52%, p < 0.001), and had worse NYHA Functional Class (median Class II vs. I, p = 0.04).Patients with diabetes had higher creatinine (96 .0]ng/mL, p < 0.001).There was no difference in LV ejection fraction and GLS between the groups.The diabetes group had a higher prevalence of LGE (76 vs. 60%, p = 0.002, driven by higher rates of ischaemic LGE (51 vs. 33%, p = 0.003).Prevalence of non-ischaemic LGE did not differ between the groups.The diabetes group had higher  After adjustment for potential confounders, including age, sex, hypertension, and coronary artery disease, diabetes remained independently associated with the presence of LGE and increased ECV (p < 0.005 for both analyses).Amongst patients without coronary artery disease, those with diabetes had a numerically higher proportion with ischaemic LGE (21%) compared to those without diabetes (13%), although this difference did not achieve statistical significance (Supplemental Table 1).

Predictors of primary outcome in all patients with heart failure
Over a median follow-up of 23 months, 74 patients experienced the primary outcome (52 hospitalisations for HF, 22 all-cause deaths; Supplemental Table 2).In univariable Cox regression analyses for the entire cohort, clinical markers associated with the primary outcome included diabetes, NYHA functional class, and systolic blood pressure.Circulating biomarkers associated with the primary outcome included NT-proBNP, hs-TnT, GDF-15, sST2, and Gal-3.CMR markers of adverse remodelling (LV mass index and end diastolic volume index), contractile function (LV ejection fraction and GLS), and myocardial fibrosis (presence of LGE, native T1 and ECV) were predictors of worse outcomes.Diabetes, NT-proBNP, and GLS remained independently associated with outcomes in the multivariable analysis (Table 2; Fig. 2).

Predictors of primary outcome in patients with heart failure and diabetes
In the diabetes group, 50 patients experienced the primary outcome (35 hospitalisations for HF, 15 all-cause deaths).In univariable analyses, systolic blood pressure, NT-proBNP, hs-TnT, GDF-15, sST2, LV ejection fraction, GLS, and native T1 were associated with the primary outcome (Table 3).GLS and sST2 remained associated with outcomes in the multivariable analysis.Patients with diabetes and worse than median GLS (GLS > -9.9%) had the worst prognosis (log-rank p < 0.001, Fig. 3).Of note, patients with diabetes and GLS better than median had similar outcomes to patients without diabetes and GLS worse than median (Fig. 2B; p = 0.70).
In a multivariable model including clinical markers and plasma biomarkers, only sST2 remained independently associated with the primary outcome, demonstrating stronger prognostic associations than NT-proBNP (Table 4).
The diabetes group was also stratified into four groups based on GLS (above or below the cohort median of -9.9%) and plasma biomarker (above or below median for NT-proBNP, hs-TnT, GDF-15, sST2, and Gal-3).Amongst patients with above median GLS, NT-proBNP and GDF-15 showed additional prognostic value (Supplemental Fig. 1).On the other hand, the combination of below median GLS and sST2 demonstrated particularly adverse prognosis (Supplemental Fig. 2).

Discussion
The main finding of this study is that GLS by CMR is a powerful, independent predictor of adverse outcomes in patients with HF and diabetes, providing incremental prognostic information beyond several circulating plasma biomarkers and CMR markers of cardiac remodelling, inflammation, and fibrosis.In the diabetes group, worse GLS (> median − 9.9%) was associated with adverse prognosis.GLS by CMR is known to have prognostic value in HF, regardless of whether the cause is ischemic or non-ischemic and whether the EF is reduced or preserved [12,13].In this study, we have extended those findings to the diabetes subgroup, who are known to have a higher risk of adverse cardiovascular outcomes.Our findings are in agreement with a prior study which showed that GLS by speckle tracking echocardiography has prognostic value in diabetes patients with dilated cardiomyopathy [26].The results from the current study raise the possibility of GLS-guided risk stratification and management in patients with HF and diabetes.Indeed, there is emerging evidence of improvement in GLS with sodium-glucose cotransporter-2 inhibitors, which may enable this in future [27,28].
GLS by CMR is well-suited for routine clinical use.The technique relies on software packages which track the endocardial and epicardial borders, and reference values are specific for the software used for analysis.Speckle tracking echocardiography has a number of limitations, including dependence on high quality 2D images, and results are often affected by foreshortening, dropout of the apical and anterolateral segments on apical views, and/or suboptimal acoustic windows [13,29,30].GLS by CMR overcomes these limitations and can be measured from routinely acquired bright-blood, steady-state free precession long-axis cine imaging.Disadvantages of GLS by CMR include its susceptibility to throughplane motion artefacts, limited temporal resolution for quantification of strain rate, and limitations in patients with contraindications to CMR [13].Furthermore, GLS by CMR may be affected by inter-vendor differences in algorithms, similar to speckle tracking echocardiography, which has resulted in difficulties establishing reference values [30].Nonetheless, GLS by CMR has shown close correlation with GLS by speckle tracking echocardiography, and has become an important component of the CMR examination alongside volumetric analysis and tissue characterisation in the assessment of HF [31].
In this study, sST2 demonstrated incremental prognostic value over other biomarkers including NT-proBNP and hs-TnT in HF patients with diabetes.sST2 is produced by cardiomyocytes and cardiac fibroblasts in response to myocardial stress, inflammation or injury [32][33][34].Other sources of sST2 include endothelial cells of the aorta and coronary arteries as well as certain immune cells [32].It acts as a decoy receptor for interleukin-33, attenuating its cardioprotective biological effects which include blunting myocardial hypertrophy and fibrosis, and inhibiting cardiomyocyte apoptosis [16,32,35].We have confirmed that sST2 levels are higher in HF patients with diabetes [36,37].Furthermore, sST2 is known to have prognostic value in acute and chronic HF, independent of plasma natriuretic peptides [38][39][40].We have provided further evidence that these findings also apply to the subgroup of HF patients with diabetes [41].Our findings therefore suggest a possible role for sST2 in risk stratification amongst HF patients with diabetes and in monitoring response to treatment [42][43][44][45][46][47].
In this study, HF patients with diabetes had worse CMR markers of myocardial injury, inflammation, and fibrosis.The association between diabetes and elevated ECV remained after adjusting for potential confounders that differed between the diabetic and non-diabetic groups including age, hypertension, coronary artery disease and presence of LGE.These findings are in agreement with prior studies, including a recent meta-analysis, which showed an association between diabetes and a higher degree of myocardial fibrosis as assessed by histology as well as ECV by CMR [9].Furthermore, diabetes remained associated with the presence of LGE in logistic regression models, even accounting for age, sex, hypertension and coronary artery disease.This was an interesting finding given that the higher rate of LGE in the diabetes group was primarily driven by higher rates of ischemic LGE.These results may reflect a higher frequency of silent myocardial infarction in diabetes patients, which has been demonstrated previously using a variety of imaging modalities, including CMR [48,49].Although our diabetes group had a numerically higher proportion with ischemic LGE in the absence of coronary artery disease, that difference did not reach statistical significance.
Both focal and diffuse myocardial fibrotic processes are known to occur in patients with diabetes independently of co-morbid conditions, including coronary atherosclerosis and hypertension.Diffuse interstitial and perivascular fibrosis are structural hallmarks of diabetic cardiomyopathy, but focal replacement fibrosis can also be seen, even in the absence of coronary artery disease [50].Several mechanisms may explain the fibrosis burden in diabetes.Hyperglycaemia is thought to upregulate the expression of profibrotic factors such as transforming growth factor beta 1 and down-regulate the activity of the matrix metalloproteinases [51].Hyperglycaemia is also known to promote the formation of advanced glycation end products (AGEs) which can cross-link collagen in the interstitium, increasing their resistance to degradation.AGEs can result in generation of reactive oxygen species and oxidative stress which further promotes a pro-fibrotic state [9].Pro-inflammatory cytokines and chemokines, as well as increased renin-angiotensin-aldosterone system activity in diabetes have also been implicated in the development of myocardial fibrosis.
Native T1 values were higher in the diabetes group and may reflect interstitial expansion due to myocardial fibrosis as well as myocardial oedema affecting the cellular and interstitial compartments [52].There is some disagreement in the literature regarding whether native T1 values are increased in diabetes.Several studies have shown an association between diabetes and increased T1 [53][54][55][56].Although a positive association was found between diabetes and native T1 in a recent meta-analysis, the result was not statistically significant [9].The authors of that study suggested that the lack of statistical significance may have resulted from limited sample size of the included studies.Indeed, the diabetes cohort in the present study was larger than any of those in the included studies.Furthermore, T1 values are known to be dependent on a variety of factors, including field strength, pulse sequence, and region of measurement within the myocardium.In the present study, T1 maps were obtained at 3T using the same MOLLI sequence for all patients, and with experienced observers performing standardized analyses that may improve reproducibility and eliminate technical cofounders [57,58].Whilst it does appear that diabetes is associated with higher T1, larger studies controlling for the variability in T1 are required to more conclusively answer this question.

Limitations
One limitation of this study is the absence of measures of glycaemic control, which are known to be prognostic in patients with HF and diabetes [59].However, GLS by speckle tracking echocardiography is known to have prognostic value independent of glycaemic control in diabetes cohorts with preserved and reduced ejection fraction, and therefore this is unlikely to have altered the prognostic associations for GLS by CMR [14,26].
Our HF cohort included patients with both ischaemic and non-ischaemic etiologies, as well as both preserved and reduced ejection fraction.Unfortunately, the study was not powered for analyses of these subgroups.Documented history of coronary artery disease did not associate with outcomes in the univariable analysis and therefore etiology of HF is unlikely to have affected the results.Finally, this was a single centre study using feature tracking software from a single vendor, limiting generalizability of the results.

Conclusions
Patients with HF and diabetes had worse CMR and plasma markers of injury, inflammation, and fibrosis, and an adverse prognosis.sST2 showed incremental prognostic value beyond NT-proBNP in HF patients with diabetes.GLS is an important and independent prognostic  marker in this group.Future studies should explore whether GLS-guided risk stratification and management can improve outcomes in this group of patients.

Fig. 2
Fig. 2 Event-free survival curves showing: (A) Adverse prognosis in patients with heart failure and diabetes; (B) Patients with diabetes and GLS worse than cohort median (-9.9%) had worst prognosis.Patients with either diabetes or GLS worse than median had similar outcomes.Abbreviations: DM: diabetes mellitus; GLS: global longitudinal strain

Fig. 3
Fig. 3 Example of GLS measurement in apical 2-chamber view.Top panel shows a patient with heart failure, diabetes and GLS ≤ -9.9%, at end systole.Bottom panel shows a patient with heart failure, diabetes and GLS > -9.9%, at end systole.Centre panel: event-free survival curves showing adverse prognosis in patients with heart failure, diabetes and GLS worse than cohort median (-9.9%).Abbreviations: DM: diabetes mellitus; GLS: global longitudinal strain

Table 2
Univariable and multivariable Cox regression models (stepwise selection) in the entire cohort for the primary outcome of allcause mortality or heart failure hospitalisation

Table 3
Univariable and multivariable Cox regression models (stepwise selection) in the diabetes group for the primary outcome of all-cause mortality or heart failure hospitalisation

Table 4
Multivariable Cox regression model (Clinical and Plasma biomarkers) in the diabetes group for the primary outcome of all-cause mortality or heart failure hospitalisation Abbreviations:; NYHA: New York Heart Association; NT-proBNP: N-terminal pro-brain natriuretic peptide; hs-TnT: high-sensitivity troponin T; GDF-15: growth differentiation factor 15; sST2: soluble ST2; Gal-3: galectin 3. § Assessed for inclusion in multivariable model but not retained